First-principles approach to Non-CoUinear Magnetism : towards Spin Dynamics 
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A description of non-collinear magnetism in tlie framework of spin-density functional theory is 
presented for the exact exchange energy functional which depends explicitly on two-component 
spinor orbitals. The equations for the effective Kohn-Sham scalar potential and magnetic field 
are derived within the optimized effective potential (OEP) framework. With the example of a 
magnetically frustrated Cr monolayer it is shown that the resulting magnetization density exhibits 
much more non-collinear structure than standard calculations. Furthermore, a time-dependent 
generalization of the non-collinear OEP method is well suited for an ab-initio description of spin 
dynamics. We also show that the magnetic moments of solids Fe, Co and Ni are well reproduced. 
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The extension of the original density functional theory 
(DFT) Hohenberg-Kohn-Sham approach to the case of 
spin polarized systems was given, under the name spin 
DFT (SDFT) more than three decades ago Q. While this 
formulation was for arbitrary directions of the magnetiza- 
tion vector field, even today most applications are based 
on a restricted collinear version. This has the advantage 
of computational simplicity: one then works with two 
separate KS equations, one yielding the spin-up orbitals 
the other the spin-down orbitals, whereas the general for- 
mulation involves Pauli spinors. Nevertheless, there ex- 
ists a wealth of non-coUinearity in nature. To give only 
a few examples it is widely seen in molecular magnets, 
exchange frustrated solids (7-Fe, spin glasses), and all 
magnets at finite temperatures. 

Crucial for practical calculations using SDFT is the ap- 
proximation made for the exchange-correlation (xc) en- 
ergy functional. The Local Spin Density Approximation 
(LSDA) and the Generalized Gradient Approximations 
(GGAs) are currently the most popular ones. These have 
been developed for collinear magnetism, and their use 
in non-collinear situations relies on the magnetization, 
m(r), and exchange correlation magnetic field, Bxc(i"), 
being made collinear in a local reference frame at each 
point in space 0]- This is only possible with purely lo- 
cal functionals like LSDA [3,3, though it has been used 
under additional approximations for gradient function- 
als as well [5]. Such approximations (that lead to lo- 
cally collinear magnetization and xc magnetic field) cause 
m(r) X Bxc(r) to vanish everywhere in space. As noted 
recently, this fact renders the adiabatic time dependent 
extension of these functionals improper for the study 



of spin dynamics, because in the absence of external mag- 
netic fields and within adiabatic approximation, the local 
torque on the spins (m(r, t) x Bxc(r, t)) vanishes This 
is a serious limitation since the dynamics of the spin de- 
gree of freedom is responsible for a number of important 
phenomena such as spin injection, the dynamics of Bloch 
walls, spin wave excitations 3, ^^'^ spin filtering, mech- 
anisms crucial for recent developments in spintronics [3]- 
The search for approximate xc functionals which depend 
on all three components of the spin magnetization m be- 
yond the form of the locally collinear LSDA has remained 
a major challenge in the description of non-collinear mag- 
netism. 

In recent years, an alternative route to the construc- 
tion of approximate xc functionals has enjoyed increasing 
interest. These involve functionals depending explicitly 
on the single-particle KS orbitals which, through the KS 
single-particle equation, are implicit functionals of the 
density [10]. Technically, one needs to employ the Opti- 
mized Effective Potential (OEP) [11] method to compute 
the local xc potential. The simplest orbital-dependent 
approximation to the xc energy is the EXact eXchange 
(EXX) functional which is the Fock exchange energy but 
evaluated with KS orbitals (i.e. orbitals coming from a 
local potential). A number of successful EXX calcula- 
tions have been reported for semiconductors 13, 13, 14 ] 
and magnetic metals [TBI ]. However, for magnetic systems 
again the collinear formalism has been employed. 

In this Letter we extend the OEP formalism for SDFT 
to non-collinear magnetic systems. Most importantly, we 
do not rely on a condition of local coUinearity and treat 
the wavefunctions as Pauli spinors for high lying and 
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Dirac spinors for deep lying (3 Ha below the Fermi level) 
electrons. Using the EXX functional, we demonstrate 
with the example of an unsupported Cr(lll) monolayer, 
that (i) the magnetization and Bxc are generally not lo- 
cally parallel in contrast to what has been assumed in 
all calculations to date and (ii) that the non-coUinearity 
is much more pronounced than found with the LSDA 
functional. Against popular belief we find that this 
non-coUinearity is not restricted to just the interstitial 
region but spreads all the way to the atom center. With 
the examples of bulk Fe, Co and Ni we further show that 
our formalism can also be effectively used for coUinear 
magnets. 

To derive the OEP equations in the general non- 
collinear case, we start with the Kohn-Sham (KS) equa- 



tion for two-component spinors which has the form of 
a Pauli equation. For non-interacting electrons moving 
in an effective scalar potential Vg and a magnetic vector 
field Bg it reads as (atomic units are used throughout) 

(^-^V^ + vs{r) + liBCT ■ B,(r) j $,(r) = £,$,(r). (1) 

This equation can be derived by minimizing the total 
energy which, in SDFT, is given as a functional of the 
density p(r) — YlT"^ magnetization 
density m(r) — /is ^"'^'^ <i>|^(r)cr$i(r). For a given ex- 
ternal scalar potential Ucxt and magnetic field Boxt this 
total energy reads 



E[p,m] = T,[p,m]+ j p{v)v^^t{v) dv + j m(r) • Bcxt(r) dr + C/[p] + i^xcb, m] 

OCC p p 

= X! ~ / P(i")^'xc(r) dv ~ \ m(r) • Bxc(r) dv - [/[p] + £'xc[p, m]. 



(2) 



where J7[p] = 1/2 JJ /3(r)p(r')/|r— r'| dr dr' is the Hartree 
energy. The xc potential and xc magnetic field are given 
by 



:(r) 



(5^xc[p, m] 
<5p(r) 



and Bxc(r) 



8E^c [p, m] 
(5m(r) 



(3) 



respectively. The exact functional form of E'xc \p^ m] is 
unknown and has to be approximated in practice. 

Assuming that the densities (p, m) are non-interacting 
(ti, B)-representable one may, equivalently, minimize the 
total-energy functional ([2]) over the effective scalar po- 
tential and magnetic field. Thus the conditions 



bE\p, m] 



^ 



and 



<5B3(r) 



= 



(4) 



must be satisfied. 

If the functional derivatives in Eq. [4] are evaluated 
for an xc functional that depends explicitly on the KS 
spinors, one obtains the natural extension of the OEP 
equations to non-coUinear magnetism. By the usage of 
spinor valued wavefunctions we can stay within a sin- 
gle global reference frame, in contrast to the case where 
functionals originally designed for coUinear magnetism 
are used in a non-coUinear context by introducing a local 
reference frame at each point in space. The most com- 
monly used orbital functional is the EXX energy given 
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El^^{{^.^\ ^ 



E 



$l(r)$,(r)<i>J(r')<i>.(r') 



dvdv' 

(5) 



where the label occ indicates that the summation runs 
only over occupied states. In the following we restrict 
ourselves to an exchange-only treatment although gener- 
alization to other orbital functionals is straightforward. 

For the energy functional Eq.© using the EXX ap- 
proximation to i?xc one obtains the following coupled in- 
tegral equations for the exchange potential and magnetic 
field 



i?.(r) 



_ 8E\pM 



(5z;s(r) 

OCC un 

EE(A 



' - ^1 



(6) 



and 



RB(r) 



8E\pM 



<SB,(r) 

OCC un / 

EE (a 



(r) 



c.c. = 0, (7) 



where Py(r) = $|(r)$j(r), my(r) = /iB$|(r)(T$j (r) 
and j runs only over the unoccupied states. The matrix 
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A is given by 

= {y^)^ - J pUr)v^ir)dr- J ml^. (r) • B.(r) dr, 

(8) 

where 

are the non-local matrix elements of the Coulomb inter- 
action between states i and j . 

To ensure that our numerical analysis be as accurate as 
possible, we use the full-potential linearized augmented 
plane wave (FP-LAPW) method ^17] implemented within 
the EXCITING code Here the single electron po- 

tential is calculated exactly without any shape approx- 
imation and the space is divided into muffin-tin (MT) 
regions, where atomic orbitals are used as a basis and 
interstitial region, where plane waves are used as a basis. 
The deep lying core states (3 Ha below the Fermi level) 
are treated as Dirac spinors and valence states as Pauli 
spinors. More importantly the magnetization density and 
xc magnetic field are both treated as unconstrained vec- 
tor fields throughout space. In our implementation of the 
OEP method the exchange fields are iteratively updated 
by subtracting the residue functions Ry and Rb from 
the exchange fields. In other words, if i is the iteration 
number then 

<(r)-«r'(r)-rK(r), 
Bi(r)=Br^(r)-rRj,(r) (10) 

is repeated until convergence is reached, with i?* and Rg 
calculated by inserting and B^~^ into Eqs. ^ and 
([7]) . r is the mixing chosen in such a manner as to achieve 
a speedy convergence. In the collinear case this method 
is similar to the one previously suggested in the Ref. [fol . 

In order to explore the impact of treating non-coUinear 
magnetism in the way outlined above we compare our 
approach with the standard LSDA functional using the 
example of an unsupported Cr (Iff) monolayer. We set 
the lattice parameter of the Cr-monolayer to that of the 
Ag (111) surface. The result is a topologically frustrated 
anti-ferromagnet, known from LSDA calculations to exist 
as a non-coUinear Neel state with the net magnetization 
direction of the three non-equivalent atoms pointing at 
120° to each other. In Fig. [1] we show the magnetiza- 
tion density and B field for both the LSDA and EXX 
functionals. Both find, as they must, the non-coUinear 
Neel state, and in fact the EXX and LSDA MT averaged 
moments are similar, being 2.60 /is and 2.0 /is, respec- 
tively. The details of the xc density and field however 
are very different with the EXX functional producing a 
lot more structure, in contrast to its fairly homogeneous 
LSDA counterpart. In the past, the LSDA results (of 
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FIG. 1: Fully non-coUinear magnetization density and B field 
obtained using the LSDA and exchange-only EXX function- 
als for an unsupported Cr-monolayer in Neel state. Arrows 
indicate the direction and information about the magnitude 
(in atomic units) is given in the colour bar. 

the kind shown in Fig. 1), which show almost no non- 
coUinearity in the MT region, led to the conclusion that 
it is sufficient to treat only the interstitial region as non- 
coUinear [3| • The present work shows that orbital func- 
tionals such as EXX are more sensitive to the atomic 
shell structure and this sensitivity also manifests itself in 
the magnetization density and exchange B field. This is 
clear from the fiower petal like structure visible in the 
magnitude of EXX density and B field. The Neel walls 
are also much narrower in the EXX case. Adding LSDA 
correlations to the EXX functional does not significantly 
change these results. A striking feature of the EXX B 
field is that, unlike its LSDA counterpart, it is not locally 
parallel to the magnetization density. 




FIG. 2: m(r) x Bx(r) for an unsupported Cr-monolayer, in 
the same plane as Fig. 1, obtained using the EXX functional. 
Arrows indicate the direction and information about the mag- 
nitude (in atomic units) is given in the colour bar. 

Another appealing property of the EXX functional 
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that could have consequences in future time-dependent 
extensions is the non- vanishing cross product of the mag- 
netization density and EXX B,, field. This is interesting 
because the equation of motion for the spin magnetiza- 
tion reads 



= 7m(r, t) X [Bxc(r, t) + Bext(r, t)] - V ■ J, 

(11) 

where Jg is the spin current and 7 the gyromagetic ratio. 
In the time-independent LSDA and conventional GGA, 
m(r) and Bxc(r) are locally collinear, as is clear from Fig. 
1, and therefore m(r) x Bxc(r) vanishes. This also holds 
true in the adiabatic approximation of time dependent 
SDFT which, by Eq. ([TT|l . implies that these functionals 
cannot properly describe the dynamics of the spin mag- 
netization. In contrast, already at the static level, for the 
EXX functional m(r) x Bx(r) does not vanish (see Fig. 
2) In fact, in the ground state of a non-coUinear ferro- 
magnet without external magnetic field, m(r) x Bxc(r) 
exactly cancels the divergence of the spin current, V • Js, 
i.e. these terms are equally important, and it is essential 
to have a proper description of m(r) x Bxc(r). These re- 
sults indicate that a time-dependent generalization of our 
method could open the way to an ab-initio description of 
spin dynamics. How well this functional really performs 
in describing the spin dynamics remains a question for 
future investigations. 

We now turn to the question of the calculation of mag- 
netic moments of collinear solids with the present formal- 
ism using the EXX functional. For the collinear mag- 
nets Fe, Co, and Ni we find moments of 2.71 (2.12 
IIb), 1.77 hb (1.71 hb), and 0.50 iib (0.55 /x^) respec- 
tively, where the LSDA results are indicated in brack- 
ets. Surprisingly, a previous OEP calculations [H, [20| 
found much larger moments of 3.40 hb, 2.25 hb and 0.68 
Hb respectively. This discrepancy may be attributed to 
the following facts: first, the previous calculations used 
the atomic sphere approximation for the scalar potential 
and the atomic moment approximation for the magne- 
tization. In our work there is no shape approximation 
for the scalar potential and the magnetization is treated 
as an unconstrained vector field. Second, and more im- 
portant, in the present work a coupled set of equations 
is solved to numerically invert the response function. 
This has the advantage of automatically including the re- 
sponse of the system to a constant magnetic field which is 
important for spin-unsaturated systems. This response 
needs additional treatment in the case where a decou- 
pled set of equations is used and the response is inverted 
in a constant-free basis, as done in all past calculations 
Til 2o| . We suspect that this is the major reason for the 
present discrepancy. 

To conclude we have presented a generalization of the 
widely used OEP equations for non-coUinear magnetic 



systems. The resulting method does not need any as- 
sumption of local coUinearity for m(r) and Bxc(r), and 
therefore extends ab-initio approaches to non-coUinear 
magnetism substantially beyond the LSDA. In particu- 
lar, a time-dependent extension of the non-coUinear OEP 
method naturally leads to a new and promising ab-initio 
approach to describe spin dynamics. 

Finally, we note that since the formalism presented 
here treats KS wavefunctions as spinors, it can be used 
in conjunction with spin-orbit coupling. In particular, in 
/-electron systems both spin-orbit coupling and the ex- 
change field are of crucial importance, where the latter is 
well known to be poorly treated by LSDA/GGA. Hence, 
the present work opens new interesting routes for future 
extensions. 
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